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Large errors in numerical weather prediction are often associated with ex- 
plosive cyclogenesis. Most studies focus on the under-forecasting error, i.e. 
cases of rapidly developing cyclones which are poorly predicted in numerical 
models. However, the over-forecasting error (i.e., to predict an explosively 
developing cyclone which does not occur in reality) is a very common er- 
ror that severely impacts the forecasting skill of all models and may also 
present economic costs if associated with operational forecasting. Unnec- 
essary precautions taken by marine activities can result in severe economic 
loss. Moreover, frequent occurrence of over-forecasting can undermine the 
reliance on operational weather forecasting. Therefore, it is important to 
understand and reduce the predictions of extreme weather associated with 
explosive cyclones which do not actually develop. 

In this study we choose a very prominent case of over-forecasting error 
in the northwestern Pacific. A 960 hPa cyclone develops in less than 24 
hour in the 5-day forecast, with a deepening rate of about 30 hPa in one 
day. The cyclone is not verified in the analyses and is thus a case of severe 
over-forecasting. 

By assimilating AIRS data, the error is largely eliminated. By following 
the propagation of the anomaly that generates the spurious cyclone, it is 
found that a small mid-tropospheric geopotential height negative anomaly 
over the northern part of the Indian subcontinent in the initial conditions, 
propagates westward, is amplified by orography, and generates a very intense 
jet streak in the subtropical jet stream, with consequent explosive cycloge- 
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nesis over the Pacific. The AIRS assimilation eliminates this anomaly that 
may have been caused by erroneous upper-air data, and represents the jet 
stream more correctly The energy associated with the jet is distributed 
over a much broader area and as a consequence a multiple, but much more 
moderate cyclogenesis is observed. 
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Abstract. Large errors in numerical weather prediction are often associated with ex- 
plosive cyclogenesis. Most studies focus on the under-forecasting error, i.e. cases of rapidly 
developing cyclones which are poorly predicted in numerical models. However, the over- 
forecasting error (i.e., the prediction of an explosively developing cyclone which does not 
occur in reality) is a very common error that severely impacts the forecasting skill of 
all models and may also present economic costs if associated with operational forecast- 
ing. Unnecessary precautions taken by marine activities can result in severe economic 
loss. Moreover, frequent occurrence of over-forecasting can undermine the reliance on op- 
erational weather forecasting. Therefore, it is important to understand and reduce the 
predictions of extreme weather associated with explosive cyclones which do not actually 
develop. 

In this study we choose a very prominent case of over-forecasting error in the north- 
western Pacific. A 960 hPa cyclone develops in less than 24 hour in the 5-day forecast, 
with a deepening rate of about 30 hPa in one day. The cyclone is not verified in the anal- 
yses and is thus a case of severe over-forecasting. 

By assimilating data from the Atmospheric Infrared Sounder (AIRS) mounted on the 
NASA’s Aqua Spacecraft, the error is largely eUniinated. By following the propagation 
of the anomaly that generates the spurious cyclone, it is found that a small mid-tropospheric 
geopotential height negative anomaly over the northern part of the Indian subcontinent 
in the initial conditions, propagates westward, is amplified by orography, and generates 
a very intense jet streak in the subtropical jet stream, with consequent explosive cyclo- 
genesis over the Pacific. The AIRS assimilation eliminates this anomaly that may have 
been caused by erroneous upper-air data, and represents the jet stream more correctly. 

The energy associated with the jet is distributed over a much broader area and as a con- 
sequence a multiple, but much more moderate cyclogenesis is observed. 


1. Introduction 

Explosive cyclo genesis, defined as a cyclone formation 
in which the deepening rate exceeds 24 hPa in 24 hours, 
normalized at the latitude of 60°, (Sanders and Gyakum, 
1980), occurs frequently in nature and is a major task 
for numerical forecasting. Over- the last 25 years, enor- 
mous improvements have been obtained in the forecasting 
of explosive cyclones. Among the several complex forc- 
ings contributing to explosive cyclogenesis (e.g Bullock 
and Gyakum, 1993), the role of latent and sensible heat 
fluxes (e.g Bo.sarfc, 1981) and the role of jet streaks (e.g. 
Uccellini et al., 1987) have received the greatest atten- 
tion from the atmospheric modeling community. Satel- 
lite information have helped substantially to improve the 
representation of surface fluxes and of the jet stream, 
the latter being related to a correct representation of the 
three-dimensional thermal structure of the atmosphere 
and to improvements in the analysis of tropopause fold- 
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mg. 

In spite of the improvements, explosive cyclogenesis 
still represents a challenge even for the current atmo- 
spheric models, and explosive cyclones can be totally 
missed in operational weather forecasts on various ranges. 
Several studies have considered the so-called 1 underfore- 
casting ’ error, in which the formation of an explosive cy- 
clone or its explosive deepening are not predicted, and 
have proposed different strategies to improve the fore- 
cast of these dangerous storms (e.g., among others, Hello 
et al., 2000; Leutbecher et al., 2002; Fourrie et al., 2003). 

However, even if the increased resolution of global 
models, the improved representation of sea surface fluxes, 
and the capability of models’ dynamical cores to develop 
rapidly growing disturbances have reduced the underfore- 
casting error, the reverse problem is being faced even 
more often: spurious, explosive cyclones develop in the 
models and are not verified in the analyses. This ‘over- 
forecasting’ error is most common when the increased 
capability of global models to represent some of the sub- 
synoptic scale aspects of explosive cyclogenesis is not 
matched by an equal improvement in the initialization. 
The more a model can develop realistically fast growing 
modes, the more important it becomes to avoid inserting 
spurious disturbances in the initial conditions. 

The overforecasting error severely affects the overall 
performance of the model, but it has not been systemat- 
ically studied. It is however a very important error, be- 
cause it undermines confidence in operational forecasting 
and is often associated with economic losses, whenever 
precautions are taken in response to a forecast of an ex- 
plosive event which does not verify. 

In this work, we select a spectacular case of overfore- 
casted cydogenesis that occurred in a 5-day forecast ini- 
tialized on 8 January 2003 over the northwestern Pacific. 
We then perform an assimilation experiment in which 
AIRS-derived thicknesses profiles are assimilated glob- 
ally, to improve the vertical temperature structure. In 
this simulation, the spurious cyclone is completely sup- 
pressed. We then compare the dynamical evolution of 
the two simulations, in order to investigate the physical 
mechanisms through which the information introduced 
by the AIRS assimilation propagates in the model. 


2. The Model 

The model used in this study is the finite- volume Gen- 
eral Circulation Model (fvGCM) developed at the NASA 
Goddard Space Flight Center (GSFC) and based on 
the finite- volume dynamical core with terrain- following 
Lagrangian control-volume discretization documented in 
Lin (2004). The finite- volume dynamical core at NASA 
GSFC is the result of more than ten years of intense ef- 
fort, in which the most fundamental steps are: 

1. Development of algorithms for transport processes 
of water vapor (Lin et al., 1994). 

2. Development of multi-dimensional “Flux-Form 
Semi-Lagrangian Transport” scheme (FFSL, Lin and 
Rood, 1996). 

3. Adaptation of the FFSL algorithm to the shallow 
water dynamical framework (Lin and Rood, 1997). 

4. Development of a simple finite- volume integra- 
tion method for computing pressure gradient in general 
terrain-followmg coordinates (Lin, 1997). 

The physics parameterizations implemented in this 
version of the NASA fvGCM axe from the the National 
Center for Atmospheric Research (NCAR) Climate Com- 
munity Model CCM version 3.6.6 (Kiehl et al., 1998). 



In tins work, the forecast model is coupled with a se- 
quential analysis system in which the atmospheric state is 
carried forward in time by the prediction model, and the 
Kalman filter analysis equation is solved at each synoptic 
time by the so-called Physical-space Statistical Analysis 
System (PSAS, Cohn et al. 1998). This system com- 
prised the fvGCM and PSAS is also called the finite- 
volume Data Assimilation System (fvDAS) and opera- 
tionally is known as the NASA Goddard Earth Observ- 
ing System (GEOS). The work described in this article is 
carried out using version 4.0.3 of the GEOS (equivalent 
to fvDAS 1.4r2), with the forecast model being run at 
a horizontal resolution of l°xl.25° and 55 layers in the 
vertical, while the analysis increments are calculated at 
a resolution of 2°x2.5° and 25 levels. 

3. The Experiment. 

3.1. A case study of overforecasting error 

The particular formulation of the finite-volume dy- 
namical core enables the model to effectively handle 
strong gradients, so that fronts and sharp thermal gra- 
dients are accurately represented and vorticity is ad- 
vected very efficiently. However, the initialization is a 
crucial aspect of the forecasting problem, and a model 
in which rapid developments axe intrinsically possible is 
often prone to large errors if the initialization is not ad- 
equate. One of the typical consequences of less-then- 
op timal initializations for a model capable of resolving 
rapid developments, is the creation of spurious rapidly 
developing cyclones. 

In this work, we select a case of severely overforecasted 
cyclogenesis that occurred over the northwestern Pacific 
in a 5-day forecast initialized on 8 January 2003. This 
simulation is defined as the “Control run” . We then per- 
form an assimilation experiment in which AIRS- derived 
thickness profiles are assimilated globally, with the aim 
of improving the overall vertical thermal structure of the 
atmosphere. We define this experiment as the “AIRS 
run”. 


3.2. Results 

In Fig. 1 the sea level pressure (sip) for the two ini- 
tializations (Control and AIRS) at 0000 UTC 8 January 
2003, together with their corresponding 120-hour fore- 
cast, axe shown. Although relatively small differences 
due to the assimilation of AIRS data are present, both 
Figs. la,c display a major, large-scale cyclone spanning 
nearly the entire northern Pacific. This system advects 
Siberian air over the Kuro-Shio and over most of the 
northwestern Pacific, and is therefore a major potential 
source of instability and secondary developments. 

In Fig. lb it can be seen that three smaller-scale cy- 
clones develop in 120 hours out of the large-scale sys- 
tem observed in Fig. la; one of these is a very intense, 
smaller-scale cyclone, with a minimum sip of 960hPa, lo- 
cated at about 42 ° N and 167 °E. The shading in Fig. lb 
represents the sip forecast minus control analyses at the 
corresponding time, 0000 UTC 13 January 2004. The 
small-scale intense cyclone does not verify in the analy- 
ses: therefore the sip forecast departure from the analyses 
is extremely large, of the order of 45 hPa over the cyclone 
center. 

Figures lc,d show that the AIRS run does not predict 
this spurious explosive cyclone. The sip difference be- 
tween the AIRS 120-hour forecast and the corresponding 
analyses at 0000 UTC 13 January is much more mod- 
est. Instead of one, very intense cyclone, the AIRS run 


produces three weak or moderate cyclones of nearly com- 
parable intensity. 

In Figs. 2 the Control and AIRS analyses at 0000 
UTC 13 January are shown, together with the NCEP 
analyses. Comparing them with Fig. 1, it can be noted 
that instead of the very intense cyclone erroneously pre- 
dicted in the Control run (Fig. lb), a very broad but 
relatively weak: cyclonic circulation spans across most of 
the northern Pacific, a remnant of the very large and in- 
tense cyclone observed in Fig. la. Within this very large 
scale cyclonic circulation, three smaller scale and rela- 
tively weak or moderate systems can be observed. It is 
a typical case in which very large-scale baroclinic condi- 
tions left by a pre-existing system can provide an envi- 
ronment favorable to intense cyclogenesis, and/or sudden 
rejuvenation, if the proper forcing occurs. However an 
erroneous prediction of any cydogenetic forcing in this 
environment can trigger rapid unstable responses in the 
model. This is what happens in the Control run, as it will 
be discussed later. On the other hand, it can be noted 
that the 120-hour AIRS forecast (Fig. Id) produces three 
smaller-scale, moderate cydone, in qualitative agreement 
with the analyses (Figs.2a,b,c). 

In Figs.2d,e, where the Control and AIRS 120-hour 
forecasts are compared with their corresponding verify- 
ing analyses at 0000 UTC 13 January over the entire 
Northern Hemisphere, it is important to observe that the 
spurious cyclone in the Control run is the largest error 
and no comparable departure can be found in the AIRS 
simulation. 

In Fig. 3 and 4 the anomaly correlation is calculated 
with respect to sip and 500 hPa geopotential height, us- 
ing the NCEP analyses for validation. It can be noted 
that AIRS improves dramatically the performance of the 
model over the Southern Hemisphere, but the impact 
is only very slightly positive over the Northern Hemi- 
sphere. This is quite remarkable and suggests that the 
standard statistical performance evaluation can occasion- 
ally be misleading. In fact in our case, smaller amplitude 
but larger- scale errors (cf. Figs. 2d and 2e) occurring 
in the AIRS run lead to a comparable global anomaly 
correlation for the two runs, even if the Control run pro- 
duces am enormous, but localized, overforecasting error. 
In this work, we focus on the individual spurious explo- 
sive cyclone in the Control forecastrand on the dynamical 
causes of its disappearance from the AIRS simulation. 

4. Dynamical considerations 

To investigate when the two simulations diverge, 
we analyze latitudinally averaged 500 hPa geopotential 
height amd sip differences with the aid of Hovmoeller di- 
agrams, on different strips of latitudes, throughout the 
northern hemisphere. The most significant result is ob- 
tained averaging the latitudes between 25 °N and 35 °N. 
In Fig. 5 the Hovmoeller diagram relative to 25 °N 
to 35 °N of 500 hPa geopotential height and sip differ- 
ence between the two simulations, AIRS and Control, is 
shown. A positive geopotential anomaly of the order of 
10m, is present in the initial conditions between longi- 
tude 80 ° E and 90 ° E (approximately over the northern 
part of the Indian subcontinent) and rapidly propagates 
eastward, doubling in amplitude between 0600 UTC 10 
and 0600 UTC 11 January. By the time the anomaly has 
reached 140° I? (which is over the Pacific Ocean), its am- 
plitude has quadrupled. Associated with this anomaly, 
two negative anomalies can be seen, moving at the same 
propagation speed but leading/lagging by approximately 
20° of longitude. The most interesting feature is a strong 



positive sip anomaly located between the maun positive 
500 hPa anomaly and the negative anomaly which is 
ahead of it. The same diagram calculated for difference 
heights shows a similar behavior between 400 and 650 
hPa approximately (not shown). 

Since we are displaying the AIRS minus Control dif- 
ference, the diagram documents that in the initial condi- 
tions, the mid-tropospheric geopotential height at about 
25 — 35° N and 80 — 90° E is lower in the Control than in 
the AIRS run. This small-scale mid-tropospheric trough 
present in the Control run amplifies rapidly between 
110 °E and 130°!?, in correspondence to the discontinu- 
ity between East Asia and the Pacific Ocean, which is 
downstream the Himalayan-Tibetan orographic forcing. 
A ridging tendency develops downstream of the moving 
trough. A sip negative anomaly indicates a cyclone in the 
Control run which should not be present, or should be 
weaker, in the AIRS run. This sip anomaly is in a tilted 
position with respect to the trough, indicating baroclinic 
growth, and rapidly propagates in phase with the trough 
and the ridge, thus benefiting from mid-tropospheric pos- 
itive vorticity advection. Since the diagram shows lati- 
tudinally averaged values between 25 °N and 35 °N, the 
magnitude of the anomalies are very significant. More- 
over, the apparent weakening at the very end of the in- 
tegration is due, as will be shown later, to a northward 
shift of the system rather than to a real weakening. 

To verify the possibility hinted by the Hovmoeller 
diagram, in Fig. 6 the 500 hPa AIRS minus Control 
height anomalies axe displayed over the areas suggested 
by Fig. 5- It can be seen that a positive (AIRS minus 
Control) anomaly is actually present over the northern 
part of the Indian subcontinent, and propagates east- 
ward throughout southern Asia The anomaly ampli- 
fies between 10 and 11 January, possibly because of the 
orographic effect of Tibet (cyclonic vorticity tends to be 
produced downstream of an orographic obstacle) . A far- 
ther deepening, not captured in the Hovmoeller diagram, 
takes place on the 12th and the 13th, when the wave 
moves northeastward, thus leaving the 25 — 35 °N aver- 
aging latitudinal band. 

According to Figs. 5 and 6, we should expect a mid- 
tropospheric shortwave developing in the Control and not 
developing (or developing less) in the AIRS simulation. 
In fact, in Fig. 7a, b we can see a small cyclonic vorticity 
maximum at 0000 UTC 9 January and at about 30—35 °N 
and 90 — 100°I£, stronger in the Control than in AIRS. 
At 0000 UTC 10 January, the vorticity maximum is at 
about 1I0°E, in agreement with Fig. 5, and at 0000 UTC 
11 January, when the strongest development starts to 
occur, the system is over the land-ocean border. At this 
time, a clear 500 hPa short wave structure is evident in 
the Control run. 

In Figs. 8 and 9 the sip, 250 hPa wind speed, 500 hPa 
geopotential height and vorticity are shown, in the Con- 
trol and AIRS simulations respectively. The longitude of 
the frame moves eastward to follow the development of 
the spurious cyclone. In the Control run, the small-scale, 
intense spurious cyclone develops in an ideal baroclinic 
way (Figs. 8c,d,e,f) in response of the mid- tropospheric 
anomaly seen before. The environment over the north- 
ern Pacific left by the decay of the very large system in 
Fig. la, is very favorable to intense baroclinic develop- 
ment, so that even a small disturbance like the one seen 
in Fig. 7e can produce a explosive cyclogenesis. It should 
be noted that the jet stream is also affected by the mid- 
tropospheric anomaly. In spite of the fact that no sub- 
stantial (AIRS minus Control) difference at 250 hPa can 
be seen in the initial conditions over the area where the 



mid- tropospheric discrepancy is found, the spurious de- 
velopment in the Control is associated with an intense jet 
streak (Fig. 8e), not present in the AIRS run (Fig. 9e). 
Moreover no evidence of the very intense 500 hPa short- 
wave can be seen in the AIRS run (cf. Figs. 8d,f with 
Figs. 9d,f) . 

5. Conclusions 

In this work we analyze a case of major overforecast- 
ing error in the NASA finite volume General Circulation 
Model. A spurious explosive cyclogenesis is produced in 
the 120-bour forecast. By comparing the model’s results 
with a simulation performed with the same model, in 
which AIRS- inferred thicknesses are assimilated globally 
and in which the spurious cyclogenesis does not occur, 
it is found that the AIRS can correct a height error that 
probably originated from incorrect upper- air observations 
over the northern part of the Indian subcontinent. This 
error causes geopotential height in the initial conditions 
to be lower of about 10-15 m and because of the partic- 
ular location (upstream of a major orographic barrier) it 
propagates as a small trough throughout Asia, is ampli- 
fied downstream by orography and, finding very favorable 
conditions over the Pacific for development, originates a 
very strong cyclone which is not observed in the analyses. 
The AIRS corrects this error in the initial conditions, and 
the corresponding simulation is thus not affected by the 
spurious development. 
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